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We  develop  the  theory  of  exchange/dipole  spin  wave  excitations  of  ferromagnetic  nanowires  of  cylindrical 
cross  section,  where  the  magnetization  is  parallel  to  the  axis  of  the  wire.  In  addition,  we  provide  the  theory  of 
the  microwave  response  of  such  structures,  for  the  case  where  the  nanowire  is  also  a  conductor.  We  present 
explicit  calculations  of  both  the  mode  structure  of  nanowires,  and  also  their  ferromagnetic  resonance  spectrum, 
with  attention  to  recent  experimental  studies.  We  compare  differences  between  the  physical  picture  appropriate 
for  the  cylinder,  with  the  well  studied  case  of  the  ferromagnetic  him. 
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I.  INTRODUCTION 

In  recent  years,  very  considerable  attention  has  been  de¬ 
voted  to  the  study  of  ultrathin  ferromagnetic  films,  and  mag¬ 
netic  multilayers  formed  from  such  films.  Such  systems  have 
unique  physical  properties,  by  virtue  of  the  fact  that  a  large 
fraction  of  the  magnetic  ions  reside  in  low  symmetry  sites  at 
the  film  surfaces,  thus  providing  strong  anisotropies  not 
found  in  bulk  magnetic  matter  constructed  from  the  same 
ions.  Also,  in  multilayers,  exchange  couplings  between  near¬ 
est  neighbor  films  provide  for  diverse  magnetic  phases,  spin 
reorientation  transitions  induced  by  very  modest  applied 
magnetic  fields,  and  other  phenomena  as  well.  Finally,  im¬ 
portant  applications  of  magnetic  multilayers  have  been  real¬ 
ized,  and  more  are  envisioned. 

Other  forms  of  magnetic  nanostructures  can  be  fabricated 
and  studied  as  well.  For  example,  Ebels  and  Wigen1  have 
created  arrays  of  very  long  ferromagnetic  nanowires  of  Ni, 
permalloy  and  Co,  with  diameters  in  the  range  of  30  to  500 
nm.  These  are  very  uniform  in  cross  section,  with  lengths  in 
the  range  of  20  microns.  They  thus  are  realizations  of 
nanowires  one  can  reasonably  view  as  infinite  in  length,  to 
excellent  approximation.  These  authors  have  carried  out  fer¬ 
romagnetic  resonance  studies  of  their  samples,  which  consist 
of  nanowire  arrays,  with  individual  entities  accurately  paral¬ 
lel  to  each  other,  but  arranged  randomly  over  a  plane. 

The  diameter  range  explored  in  these  studies  is  such  that 
when  one  considers  the  spin  wave  modes  of  an  isolated 
nanowire,  it  is  necessary  to  include  both  exchange  and  dipo¬ 
lar  contributions  to  their  excitation  energy.  The  present  paper 
presents  the  theory  of  such  dipole/exchange  spin  wave 
modes  of  isolated  ferromagnetic  nanowires  of  circular  cross 
section,  for  the  case  where  the  magnetization  is  parallel  to 
the  symmetry  axis  of  the  wire.  In  addition,  we  provide  a 
description  of  the  ferromagnetic  resonance  absorption  spec¬ 
trum,  for  the  case  where,  as  in  the  samples  discussed  above, 
the  materials  are  metallic  in  nature.  It  is  thus  important  to 
take  into  account  the  influence  of  the  finite  skin  depth,  par¬ 
ticularly  for  the  larger  diameters  explored  in  these  particular 
studies. 

There  is  a  long  history  of  the  study  of  the  magnetostatic 


spin  wave  modes  of  uniformly  magnetized  ferromagnets  of 
various  shapes,2  with  explicit  attention  to  cylindrical 
samples.3  However,  little  attention  has  been  devoted  to  the 
inclusion  of  exchange,  since  these  early  studies  were  moti¬ 
vated  by  applications  to  ferrite  samples  of  rather  large  di¬ 
mensions.  The  structure  of  the  spin  wave  modes  accessible 
to  microwave  excitation  in  such  samples  are  influenced  little 
by  exchange.  In  the  particular  case  of  samples  of  cylindrical 
shape,  one  study  of  the  exchange  dipole  modes  has 
appeared;4  so  far  as  we  know,  the  text  is  not  available  in 
English,  so  its  content  is  of  limited  accessibility. 

There  are  two  interesting  implications  of  the  mode  struc¬ 
ture  of  the  ferromagnetic  nanowires.  First,  as  we  shall  see 
below,  at  long  wavelengths,  the  dispersion  is  controlled  by 
the  long  ranged  dipolar  fields.  While  details  are  different, 
very  much  as  in  ultrathin  films,5'6  the  dipolar  contribution 
can  produce  downward  curvature  in  the  dispersion  relation  at 
long  wavelengths.  At  larger  wave  vectors,  the  curvature  is 
always  positive,  as  a  consequence  of  the  presence  of  ex¬ 
change.  Thus,  one  realizes  short  wavelength  spin-wave 
modes  degenerate  in  frequency  with  the  uniform  mode  of  the 
cylinder  excited  in  ferromagnetic  resonance.  If  the  cylinder 
surfaces  are  not  perfectly  smooth,  the  two  magnon  dipolar 
mechanism  operative  in  ultrathin  films6'7  should  then  also  be 
“active”  in  the  ferromagnetic  nanowires.  This  will  lead  to 
extrinsic  contributions  to  the  ferromagnetic  resonance  line- 
width  similar  to  those  demonstrated  to  be  substantial  in  ul¬ 
trathin  ferromagnetic  films.7 

Also,  when  the  ferromagnetic  resonance  mode  of  the 
nanowire  is  excited,  the  precession  of  the  magnetization 
leads  to  large  magnetic  fields  outside  the  sample,  in  contrast 
to  the  case  of  thin  films,  where  the  field  is  confined  to  the 
film  interior.  There  should  then  be  strong  interactions  be¬ 
tween  nanowires,  in  dense  arrays  such  as  those  explored  by 
Ebels  and  Wigen.1  The  consequence  of  these  interactions  is 
not  explored  here,  though  our  formal  structure  can  be  used 
for  this  purpose. 

In  the  discussion  above,  we  have  discussed  on  the  experi¬ 
mental  studies  of  Ebels  and  Wigen. 1  It  should  be  remarked 
that  Demokritov  and  Hillebrands8  have  carried  out  extensive 
studies  of  nanowire  arrays.  Their  nanowires  have  rectangular 
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cross  section,  however.  It  would  be  highly  desirable  to  de¬ 
velop  a  theory  of  the  dipole/exchange  modes  of  such  entities. 
We  shall  address  such  questions  in  future  studies. 

The  outline  of  this  paper  is  as  follows.  In  Sec.  II,  we 
present  the  theory  of  the  exchange/dipole  spin  wave  modes, 
and  in  Sec.  HI  we  provide  numerical  studies  of  the  spin  wave 
dispersion  and  ferromagnetic  resonance  response.  We  also 
discuss  differences  between  the  mode  structure  and  response 
characteristics  of  the  nanowire,  and  the  well  known  case  of 
the  thin  film.  In  regard  to  the  microwave  response,  for  the 
case  where  the  sample  is  made  of  conducting  material,  there 
are  very  substantial  differences  between  the  two  cases.  Sec¬ 
tion  IV  is  devoted  to  concluding  remarks. 

II.  THEORETICAL  DISCUSSION 

In  this  section,  we  first  present  the  theory  of  the  spin- 
wave  excitations  of  a  long  ferromagnetic  nanowire  with  cir¬ 
cular  cross  section,  and  magnetization  parallel  to  its  axis. 
Then  we  turn  to  a  description  of  its  response  to  a  spatially 
uniform  microwave  field,  appropriate  to  a  ferromagnetic 
resonance  experiment.  Throughout  this  section,  the  magneti¬ 
zation  is  assumed  parallel  to  the  z  axis,  which  is  also  parallel 
to  the  symmetry  axis  of  the  nano  wire. 

As  mentioned  in  Sec.  I,  we  develop  here  the  theory  of  the 
spin  waves  of  the  nanowire,  with  emphasis  on  the  regime 
where  both  exchange  and  dipolar  couplings  between  the 
spins  are  comparable  in  magnitude.  However,  before  we  turn 
to  this  general  analysis,  it  will  be  useful  for  what  follows  to 
have  in  hand  a  summary  of  the  theory  of  dipolar  spin  waves, 
in  the  absence  of  exchange.  This  will  serve  also  to  introduce 
notation  used  throughout  the  paper. 

In  the  linearized  version  of  spin  wave  theory,  when  a 
mode  is  excited  the  magnetization  of  the  sample  is  given  by 

M(r,f)  =  £M  s+ xm  X(r)exp{  —  iilt)  +  ym3,(r)exp(  —  /Of), 

(1) 

where  for  any  nanowire  of  uniform  cross  section, 
m a(x ,y ,z)  =  m a(x ,y)exp(ikz),  where  k  is  the  wave  vector  of 
the  mode,  in  the  direction  parallel  to  the  axis  of  symmetry. 
The  precession  of  the  magnetization  generates  a  magnetic 
field  of  dipolar  origin,  with  frequency  17.  We  call  this  Ir. 
The  dependence  of  the  dipolar  field  on  z  is  the  same  as  that 
given  above.  In  the  absence  of  exchange,  the  components  of 
the  dipolar  field  and  those  of  the  transverse  magnetization 
are  linked  by  the  constitutive  relations,  suppressing  explicit 
reference  to  time  dependence  for  brevity. 


hd(r)=-V<DM(r),  (3) 

where  again  $M(r)  =  <f>M(x,y)exp(fe).  Inside  the  nanowire, 
the  induction  b  =  h'/+4Trm,  while  outside  b  =  h'/.  The  mag¬ 
netostatic  theory  of  spin  waves,  with  exchange  ignored,  fol¬ 
lows  by  exploring  the  solutions  of  V  ■  b  =  0  everywhere,  with 
the  solution  subject  to  the  condition  that  the  normal  (radial, 
in  this  case)  component  of  b  be  continuous  across  the  sur¬ 
face,  while  tangential  (azimuthal,  z)  components  of  hf/  are 
continuous.  Inside  the  cylinder,  the  magnetic  potential  satis¬ 
fies  Walker’s  equation2 

/  3‘ ^  \ 

[  1  +  4  TTXi  ( ft )]  |  ^2  +  ^2)  ®M(x,y )  -  k2®M(x,y )  =  0 

(4) 

while  outside,  one  sets  ^(17)  to  zero.  We  let  pl(il)  =  1 
+  4ttxi(Q)  and  p,2(17)  =  4Tr^2(f7)  in  what  follows. 

For  the  problem  of  interest,  we  use  cylindrical  coordi¬ 
nates,  and  seek  solutions  for  which  the  magnetic  potential 
has  the  form  <f>M(p,  </>)=/m(p)exp(»H(/>).  We  shall  focus  our 
attention  here  on  the  frequency  regime  where  p.^17)  is  posi¬ 
tive.  Then  inside  the  nanowire,  one  has  fm(p) 
=AIm[kp//xl(i 7)1/2]  while  outside  one  has  />(p) 
=  BKnl(kp),  where  /,„  and  Km  are  the  modified  Bessel  func¬ 
tions. 

Application  of  the  boundary  conditions  described  above 
at  the  surface  of  the  nanowire,  where  p  =  R ,  leads  one  to  the 
implicit  dispersion  relation  from  which  the  frequency  of  the 
magnetostatic  modes  are  found: 

m\UkRlnT)\  mfi2_K'JkR ) 

/t‘  [iJkRI^2)]  kR  KJkR )■ 

In  this  expression,  I'm  and  K'm  are  derivatives  of  the  modified 
Bessel  functions  with  respect  to  their  argument. 

For  what  follows,  our  interest  will  reside  in  the  behavior 
of  the  magnetostatic  modes  at  very  long  wavelengths,  in  the 
regime  kR  <§  1 .  Through  use  of  the  appropriate  series  expan¬ 
sions  for  the  modified  Bessel  functions  in  Eq.  (5),  one  may 
obtain  analytic  expressions  for  the  frequencies  of  the  various 
modes.  We  denote,  for  a  given  choice  of  the  wave  vector  k 
and  the  azimuthal  quantum  number  m  the  frequency  of  the 
mode  by  17 m(k).  One  then  has 


1 17  j (k) |  =  —  Cll(k)  =  H0+2TTMs  —  n tM s(kR)2  In 


kR 


(6a) 


mx(x,y)  =  X\Whdx(.x,y)  +  ix2{tt)hd(x,y)  (2a) 

and 

my{x,y)  =  Xi{&)hdy{x,y)- ixi{&)hdx{x,y),  (2b) 

where  for  a  ferromagnet,  Xi ( 11)  =  !7W!7W /( lip  — 11 2)  and 
^2(0 )  =  — 172).  If  y  is  the  gyromagnetic  ratio 

and  H0  an  applied  dc  magnetic  field  parallel  to  the  magneti¬ 
zation,  we  have  introduced  17 H=yH0  and  11M=  yM s . 
Within  the  magnetostatic  approximation,  the  dipole  field 
may  be  expressed  as  the  gradient  of  the  magnetic  potential 


while  for  m  >  1 

1  1  ttMs 

\ttm{k)\=-(Dm{k)  =  H{)  +  2iTMs-  — 2— yy  (£/?)“  +  •••  . 

(6b) 

These  results  require  some  comment.  First  of  all,  when  the 
azimuthal  quantum  number  is  positive,  as  assumed  in  deriv¬ 
ing  Eqs.  (6),  the  frequencies  are  all  negative,  as  one  sees 
from  these  results.  These  modes  all  describe  a  circulation  of 
energy  in  the  counterclockwise  sense  around  the  cylinder,  as 
one  looks  down  the  magnetization.  If  we  had  chosen  the 
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azimuthal  quantum  number  negative,  the  frequencies  would 
have  all  been  positive;  the  frequency  is  an  odd  function  of  m 
as  one  can  appreciate  from  the  structure  of  Eq.  (5).  Thus  the 
entire  spectrum  of  magnetostatic  modes  describes  a  circula¬ 
tion  of  energy  around  the  magnetization,  with  the  counter¬ 
clockwise  sense. 

As  the  wave  vector  k— >0,  all  the  magnetostatic  modes 
approach  the  frequency  H0  +  2ttM s  .  This  is  the  ferromag¬ 
netic  resonance  frequency  of  the  cylindrical  nanowire.  When 
this  uniform  ferromagnetic  resonance  mode  is  excited,  the 
transverse  magnetization  engages  in  circular  precession 
about  the  static  magnetization,  and  generates  an  internal, 
spatially  uniform  demagnetizing  field  of  strength  2  itM s  as  it 
does  so.  As  we  move  off  to  finite  wave  vectors,  the  magnetic 
dipole  interactions  produce  a  negative  dispersion  initially,  so 
the  frequency  of  the  spin  wave  modes  drops  below  that  of 
the  ferromagnetic  resonance  mode.  We  will  see  below  that 
when  exchange  is  added,  a  positive  contribution  to  the  dis¬ 
persion  is  produced  for  many  of  the  modes,  so  the  minimum 
frequency  can  lie  away  from  k  =  0. 

A  complete  description  of  the  spin-wave  modes  of  the 
cylinder  requires  us  to  include  the  influence  of  exchange.  We 
now  turn  to  this  question. 


A.  Spin  wave  dispersion  in  the  presence  of  exchange 


larized  variables  m  +  _  =  m x±  imy ,  and  similarly  for  the  di¬ 
pole  field  components.  Then  Eqs.  (7)  are  replaced  by 


(D.  +  H0  —  DV2)m  +  =  M  shd+  ,  (8a) 

(il  —  H0  +  DV2)m-=  —Mshd_  ,  (8b) 


7  d 

d  \ 

d 

d  \ 

— 

'  —  \m  ,  + 

- h 

=  0. 

\  dx 

dy) 

dx 

dy) 

(8c) 


In  what  follows,  we  shall  utilize  the  identities 


d  d  \ 
dx  1  dy) 


Vi  $m. 


(9) 


where  V2  =  (d2/ dx2  +  d2l  dy2) .  We  next  introduce  the  auxil¬ 
iary  quantities 


/+ 


d  d  \ 
dx~l  dyjm  + 


(10a) 


and 


Id  d\ 

<10b) 

Then  one  may  rewrite  Eq.  (8c)  to  read 


In  the  macroscopic  description  of  the  spin  excitations,  the 
magnetic  field  h,/  inside  and  outside  the  cylinder  is  still  de¬ 
scribed  through  introduction  of  a  magnetic  potential  as  in  Eq. 
(3),  and  we  still  require  that  Vb  =  0  everywhere,  just  as 
before.  However,  we  no  longer  can  utilize  Eqs.  (2)  to  relate 
the  transverse  magnetization  components  to  the  dipolar  field. 
We  must  resort  to  the  Landau-Lifshitz  equation  of  motion 
instead,  and  in  this  one  incorporates  the  exchange  field  a 
given  spin  experiences  from  its  neighbors.  This  may  be  done 
by  replacing  the  spatially  uniform  dc  magnetic  field  by  the 
effective  field  z(H0  —  DW2),  where  D  is  the  exchange  stiff¬ 
ness.  This  leads  us  to  the  following  relations: 


i£lmx=(H0  —  DV2)my+Ms — - — ,  (7a) 

—  i£lmy  =  (H0— D\2)mx—M  s — - —  (7b) 


while  inside  the  material,  the  V  •  b  =  0  condition  becomes 


dm,.  dm,, 
V2®m- 4771-^+-^ 
ox  dy 


=  0. 


(7c) 


Outside  the  material,  the  magnetic  potential  satisfies 
Laplace’s  equation.  Of  course,  when  the  exchange  stiffness 
is  set  to  zero,  Eqs.  (7a)  and  (7b)  produce  results  equivalent  to 
those  in  Eq.  (2). 

Some  manipulation  is  required  to  cast  Eqs.  (7)  into  a  form 
where  one  may  extract  the  structure  of  the  solution,  for  the 
geometry  of  present  interest.  We  first  rewrite  the  equations, 
to  express  them  in  terms  of  the  right  and  left  circularly  po- 


/++/-  =  ^V2cbM,  (11) 

Z  77 

while  Eqs.  (8a)  and  (8b)  may  be  rearranged  to  state 

n(/++/-)  +  (H0-DV2)(/+-/_)  =  0  (12a) 


«( f+~f-)  +  (H0-DV2)(f++f_)=-2MsVl^M . 

(12b) 

Upon  combining  Eq.  (12a)  with  Eq.  (12b),  one  may  relate 
(/+  — /_)  to  the  potential  (I>  M  .  When  this  statement  and  Eq. 
(11)  are  substituted  into  Eq.  (21a),  we  obtain  a  homogeneous 
equation  that  must  be  satisfied  by  the  magnetic  potential 

[(DV2-H0){DV2-B0)]V2®m 

,  d2 

+  47rMs(DV2-//0)-TcbM  =  o.  (13) 
dz 

In  Eq.  (13),  we  introduce  B0  =  H0  +  4ttM s  . 

One  finds  that  Eq.  (13)  admits  solutions  of  the  form 


<&M(p,(/),z)  =  Jm(Kp)exp(im(/)+ikz).  (14) 

If  Eq.  (14)  is  inserted  into  Eq.  (13),  and  one  notes  that 


1 


d  d  I 
P~  +  1 


K2  — 


J  (X 


(15) 


p  dp  dp  \  p" 

then  for  Eq.  (14)  to  be  a  solution  of  Eq.  (13),  k  must  be  a 
root  of 


D2(k2  +  k2)2  +  D(H0  +  B  0)(k2  +  k2)2 

+  (H0BQ-Cl2-4TTMsDk2){K2  +  k2)-4iTMsH0k2  =  0. 

(16) 
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Since  Eq.  (16)  is  a  cubic  equation  in  k1,  for  each  choice  of  m 
and  k,  there  are  three  linearly  independent  solutions  of  Eq. 
(14).  Thus,  we  write  the  magnetic  potential  inside  the  mate¬ 
rial  in  the  form 


3 

0>M(p,^,z)  =  X  AiJm{Kip)exp(im(f>+ikz).  (17) 

;  =  i 

Outside  the  cylinder,  the  magnetic  potential  satisfies 
Laplace’s  equation  precisely  as  in  the  magnetostatic  theory, 
so  once  again  for  p>R  we  have 

<&M(p,<l>,z)=BKm(kp)exp(im(f>+ikz).  (18) 

To  find  the  dispersion  relation  of  the  spin  waves,  through 
application  of  the  boundary  conditions,  we  shall  obtain  four 
homogeneous  equations  for  the  four  coefficients  in  the  mag¬ 
netic  potential,  as  displayed  in  Eqs.  (17)  and  (18).  Two  of  the 
boundary  conditions  are  continuity  of  the  magnetic  potential 
(this  assures  continuity  of  tangential  components  of  hr/)  and 
continuity  of  bp .  We  shall  have  two  other  boundary  condi¬ 
tions  on  the  transverse  components  of  magnetization,  stated 
below.  To  apply  the  boundary  conditions,  within  the  material 
we  need  explicit  expressions  for  the  radial  and  azimuthal 
components  of  the  transverse  magnetization.  To  obtain  these, 
we  expand  these  components  in  a  Bessel  function  series  and 
use  Eqs.  (8a)  and  (8b),  noting  Eq.  (15).  The  following  op¬ 
erator  identities  are  useful: 


and 


d  d  .Jdid 
dx  dy  \dp  p  d(f> 

d  d  .Jdid 

- i  —  =e-“t’\  —+ - 7 

dx  dy  \dp  p  dcf> 


(19a) 


(19b) 


One  finds  the  following  expressions  for  the  two  components 
of  magnetization: 


m+(p,<f>,z) 


Ms  l  y  KiAiJm+1(Kip)  \ 
4  IT  \  ft  [D(J  +  k2)+HQ+kl]j 


Xexp  [i(m+l)<f>+ikz]  (20a) 


and 


m_(p,  4>,z) 


M_s_  I  y  KjAjJ m-\(Kip)  \ 

4  rr  ( f=i  [D(Kr  +  k2)  +  HQ-(l]j 


Xexp[/(f»—  \)cf>+ikz]  (20b) 

The  two  boundary  conditions  which  supplement  those  of 
magnetostatics  refer  to  the  behavior  of  the  magnetization  at 
the  cylinder  surface.  If  one  considers  an  idealized  Heisen¬ 
berg  ferromagnet,  and  examines  the  behavior  of  the  spin 
wave  eigenfunction  at  the  sample  surface,  in  the  long  wave 
length  limit,  it  is  well  know  that  one  must  require  that  the 
normal  derivative  of  the  transverse  components  of  magneti¬ 
zation  vanish  there.  We  wish  to  consider  the  possibility  that 
surface  anisotropy  is  present  on  the  cylinder  surface;  its  pres¬ 
ence  influences  the  boundary  condition.  We  assume  that  the 


surface  anisotropy  contributes  a  term  to  the  surface  energy  of 
our  cylinder  in  the  form  — Ks(mp/M s )2,  where  the  units  of 
Ks  is  ergs/cm2.  When  Ks  is  positive,  the  normal  to  the  sur¬ 
face  is  an  easy  axis,  and  when  it  is  negative,  the  normal  to 
the  surface  is  a  hard  axis.  In  terms  of  the  standardly  defined 
exchange  stiffness  used  often  in  the  literature  on  ferromag¬ 
netism!,  then  the  boundary  conditions  for  the  transverse 
magnetization  become 


dm4 
.  dp 


dm f 
dp  I 

p  '  p=R 


=  0, 

(21a) 

p  =  R 

2Ks\  , 

(21b) 

The  exchange  stiffness  D  in  our  formulas  has  units  of 
Ga-cm2.  It  is  proportional  to  A,  and  in  fact  D  =  2AI  p>BM  s , 
where  p.B  is  the  Bohr  magneton. 

This  completes  our  discussion  of  the  formalism  for  the 
description  of  the  exchange/dipole  spin  wave  excitations  of  a 
ferromagnetic  nanowire  of  cylindrical  cross  section.  We  de¬ 
scribe  numerical  calculations  based  on  this  description  in 
Sec.  III.  We  turn  next  to  the  theory  of  the  microwave  re¬ 
sponse  of  a  conducting  ferromagnetic  nanowire,  before  we 
present  the  numerical  results. 


B.  The  microwave  response  of  an  isolated  ferromagnetic 
nanowire 

In  this  section,  we  present  the  theory  of  the  microwave 
response  of  the  nanowire  considered  above,  with  exchange 
and  the  presence  of  surface  anisotropy  included.  The  geom¬ 
etry  is  the  same  as  that  considered  in  the  previous  section. 
We  have  a  nanowire,  whose  saturation  magnetization  is  di¬ 
rected  parallel  to  its  symmetry  axis.  It  is  the  case  as  well  that 
the  samples  employed  in  the  experiments  which  motivated 
this  study  are  metallic  in  nature,  so  we  wish  to  include  the 
influence  of  the  conductivity  of  the  wire  on  its  response.  In 
essence,  the  microwave  fields  to  which  the  wire  is  exposed 
create  eddy  currents  which  produce  a  finite  skin  depth.  This 
influences  the  spatial  profile  of  the  exciting  field  within  the 
sample,  though  the  discussion  presented  in  Sec.  Ill  will  show 
that  the  influence  of  the  skin  depth  is  much  more  modest  for 
the  cylindrical  geometry,  compared  to  the  ultrathin  film. 

In  what  follows,  it  is  assumed  that  the  microwave  mag¬ 
netic  field  of  interest  is  spatially  uniform  far  from  the  wire, 
and  is  parallel  to  the  x  direction.  All  quantities  will  be  as¬ 
sumed  independent  of  z  in  what  follows.  Such  a  microwave 
field  excites  the  magnetization  of  the  sample,  of  course,  so 
the  field  is  inhomogeneous  near  the  wire,  but  always  in  the 
plane,  and  independent  of  z.  By  Faraday’s  Law,  the  time 
varying  magnetic  field  generates  an  electric  field,  parallel  to 
z.  It  is  this  electric  field  which  is  responsible  for  the  eddy 
currents  in  the  conducting  material  that  lead  to  the  skin  ef¬ 
fect. 

We  may  describe  these  fields  by  introducing  the  vector 
potential 


A(x,y,t)=zA(x,y)exp(  —  iklt).  (22) 
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Then,  upon  dropping  explicit  reference  to  the  time  depen¬ 
dence  of  various  quantities  once  again,  the  electric  field  is 
given  by 

a 

e(x,y)  =  i  —zA(x,y).  (23) 

c 

The  magnetic  field  h  is  linked  to  the  electric  field  via 
4  tt  4  7 tct  2  i 

V  Xh= — j= - ze=~zrzA,  (24) 

c  c  Sq 


I  dm..  dm..\  .  . 

i^  —  +^L\=[Bo-D^2JF  +  Ms^lA.  (29) 

The  combination  of  Eqs.  (25)  and  (27)  admit  solutions  of  the 
form 


A  =aJm{Kp)exp(im<fi) 

(30a) 

F  =  bJm(Kp)exp(im(f)). 

(30b) 

where  cr  is  the  conductivity  of  the  material.  We  have  intro¬ 
duced  the  classical  skin  depth  80  =  cI(2,tt(jCI)112.  In  Eq. 
(24),  we  ignore  the  displacement  current  term  in  the  Max¬ 
well  equation.  If  this  were  to  be  included,  its  influence  can  be 
absorbed  into  a  correction  to  the  formula  for  the  skin  depth. 
The  correction  is  of  no  quantitative  importance  at  microwave 
frequencies,  for  typical  metals. 

The  magnetic  induction  b  =  V  X  A,  and  of  course  h  =  b 
—  47rm.  When  these  statements  are  combined  with  Eq.  (24), 
and  it  is  noted  that  all  quantities  are  independent  of  z,  we  find 


7  2*1 

dmx  dmy 

Vl+T2 

"0. 

£ 

II 

dy  dx 

where  the  operator  V2  is  introduced  just  after  Eq.  (9).  Addi¬ 
tional  relations  between  the  magnetization  and  the  vector 
potential  follow  from  the  Landau  Lifshitz  equations.  We  use 
analogs  of  Eq.  (7),  where  now  the  field  h=VxA— 47rm, 
rather  than  as  earlier,  and  we  also  add  to  the  right 

hand  side  the  damping  term  ( G/Mf)(M.XdM/dt ),  with  G 
the  Gilbert  damping  factor.  We  then  find,  with  all  quantities 
once  again  expressed  in  magnetic  field  units 


iClmx= 


[■ B0-DVl]my  +  Ms 


dA 

dx 


(26a) 


and 


iflmy=  ~\B0—DV±]mx+Ms— — .  (26b) 

We  have  introduced  B0  =  H0  +  4ttM  s+  ig£l,  where  g 
=  G/yMs.  After  a  bit  of  manipulation,  from  Eq.  (26),  we 
may  obtain  the  statement  that 

[(B0-DVi)2-n2](^-  -Ms(B0-DVl)VlA 

=  0,  (27) 

which  when  combined  with  Eq.  (25)  provides  us  with  two 
equations  that  link  the  vector  potential  with  the  quantity 


F(x,y) 


dmx  dmy 
dy  dx 


(28) 


For  constructing  explicit  expressions  for  the  components  of 
transverse  magnetization,  it  is  useful  to  note  the  relation 


When  these  forms  are  substituted  into  the  two  equations,  we 
find  that  k  satisfies 


D2k6  + 


2  iD2 

2{Hq  +  2ttMs)D  -=j— 

"o 


k4  + 


~  _  4  iDB0 

BqHq  “2 
°0 


2 i 

-T[JB2-n2]=o. 


(31) 


Here  we  have  H0  =  H0  +  igD,.  Notice  in  the  limit  that  the 
skin  depth  is  allowed  to  become  infinite  and  also  when 
damping  is  ignored,  Eq.  (31)  reduces  to  Eq.  (13)  if  there  we 
set  k  equal  to  zero. 

There  are  three  independent  solutions  of  Eq.  (31)  for  k. 
Thus,  the  most  general  solution  for  the  vector  potential  in  the 
medium,  with  the  azimuthal  variation  as  given  in  Eq.  (30a), 
is 


3 

A{p,<f>)  =  2  &iJm(Kip)exp(im<f>).  (32) 

i  =  i 


For  the  purpose  of  applying  boundary  conditions  at  the  sur¬ 
face  of  the  cylinder,  we  require  explicit  expressions  for  the 
magnetization,  and  for  the  components  of  magnetic  field 
generated  by  the  motion  of  the  magnetization.  It  is  most 
convenient  for  these  to  be  expressed  in  cylindrical  compo¬ 
nents.  After  some  algebra  which  employs  standard  Bessel 
function  identities,  we  find 


n(p,4>)  =  i  —  Zj  UtKt 


F m+  1  (  ^/P) 


2  i=  1  '|  b0+dk2+ n 


J m—  1  (  KiP)  .. 

H - 1  exp  (imcp), 

B0  +  DK2-fl  J 


(33a) 


,  .  vi  Jm+\  (KiP) 

*^(P,<P)=  —  Zj  CliKA  - ; - 

2  '■='  {B0  +  Dk2  +  CI 


F m  —  1  (  ^iP) 

B  o  +  D  kJ  — 17 

3 


exp  (33b) 


l  /  /  L  %  ^ 

b  {p.(/))=  Zj  aiJ  mi  KiP)  exP(  * m  (f>) ,  (34a) 

P  1=1 
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b^pA) 


A 2  aiKi{Jm  +  l(Kip) 
4  i=  l 


III.  NUMERICAL  STUDIES  OF  THE  EXCHANGE/DIPOLE 
SPIN  WAVE  MODES  AND  THE  MICROWAVE 
RESPONSE  OF  FERROMAGNETIC  NANOWIRES 


(35a) 

and  the  electric  field  is  given  in  Eq.  (23).  We  now  need  to 
match  the  interior  solutions  just  described  to  the  fields  out¬ 
side  the  cylindrical  nanowire.  It  is  a  simple  matter  to  gener¬ 
ate  expressions  for  the  magnetic  and  electric  field  in  the 
vacuum  outside,  following  a  simplified  version  of  the 
method  just  given.  One  may  generate  an  expression  for  the 
vector  potential  in  spherical  coordinates  through  the  wave 
equation  in  vacuum,  then  use  this  to  find  the  components  of 
the  various  fields.  We  do  this  assuming  the  limit  ilR/c<  1 
applies  for  the  nanostructures  of  interest;  in  the  near  vicinity 
of  the  nanowire,  the  incident  microwave  field  is  viewed  as 
spatially  uniform,  with  magnitude  hQ,  and  also  parallel  to  the 
x  axis.  We  have  outside  the  cylinder  the  incident  field, 
supplemented  by  that  generated  by  the  motion  of  the  magne¬ 
tization  in  the  nanowire;  the  latter  will  vary  with  the  azi¬ 
muthal  angle  as  exp (±i(f>)  in  the  limit  the  incident  field  is 
approximated  as  spatially  uniform  outside  the  wire.  Then  for 
the  fields  outside  we  find 


h>(p,cl>)  =  hocos(<f>)+i 


b  + 

—  exp  {icf>)  +  i 


b  _ 

—  exp  (-!», 

(36a) 


>/  b+  b- 

h4,(pA)  =  h0sm((t>)+  —rexp(i(!>)+  -^Texp(-i(f>) 

(36b) 

and  for  the  magnitude  of  the  electric  field,  we  have 

n [  b+  b_  1 

e(p, <j>)  =  i  —  h0p  sin(0)H - exp(/ <^>) H - exp(  —  icf>)  . 

c  l  P  P  ) 

(36c) 

Our  next  task  is  to  apply  the  boundary  conditions.  The  solu¬ 
tions  inside  the  medium  are  formed  by  superimposing  the 
m  =  +  1  and  m  =  —  1  forms  given  above.  Thus,  we  have  six 
unknown  amplitudes  associated  with  the  solution  in  the  inte¬ 
rior  of  the  wire,  and  as  one  sees  from  Eq.  (36)  there  are  two 
more  associated  with  the  fields  outside,  for  a  total  of  8.  The 
boundary  conditions  on  the  fields  are  of  course  that  the  mag¬ 
nitude  of  the  (tangential)  electric  field  be  conserved  across 
the  boundary,  along  with  continuity  of  bp  and  h ,,, .  The  con¬ 
servation  of  bp  and  tangential  electric  field  yield  identical 
constraints,  so  from  the  field  conservation  conditions,  we 
have  two  independent  statements.  We  also  apply  the  bound¬ 
ary  conditions  on  the  magnetization  given  in  Eq.  (21).  When 
these  four  statements  are  broken  down,  with  components 
proportional  to  and  components  proportional  to  exp( +  /</>) 
separated,  we  obtain  eight  inhomogeneous  equations  from 
which  the  unknown  amplitudes  may  be  expressed  in  terms  of 
h0.  We  will  not  quote  the  explicit  form  of  these  statements 
here,  since  they  may  be  derived  readily  from  the  information 
given  above.  As  remarked  earlier,  in  Sec.  Ill,  we  present  a 
series  of  numerical  studies  based  on  the  formalism  just 
given. 


A.  The  nature  of  the  spin  wave  modes 

First,  we  begin  with  some  general  remarks.  An  interesting 
question  is  the  initial  behavior  of  the  spin  wave  dispersion 
curves  in  the  limit  of  small  wave  vector  k.  If  initially,  there  is 
downward  curvature,  then  the  minimum  frequency  of  a  given 
branch  will  be  at  a  finite,  nonzero  wave  vector,  since  we  can 
expect  exchange  to  dominate  at  large  values  of  the  wave 
vector.  Conversely,  if  the  initial  curvature  is  positive,  then 
we  may  expect  the  minimum  to  reside  at  zero  wave  vector. 

To  examine  this  behavior,  as  one  sees  from  Eqs.  (6),  we 
need  to  consider  the  case  m=  1  separately  from  the  case 
where  m>  1.  As  we  have  seen  from  Sec.  II B,  the  zero  wave 
vector  modes  with  m  =  1  are  the  modes  excited  by  a  spatially 
uniform  microwave  field. 

We  may  expect  exchange  to  add  a  term  to  the  dispersion 
relation  proportional  to  Dk2,  of  course.  If  we  then  examine 
Eq.  (6a),  we  see  that  by  virtue  of  the  logarithmic  term,  in  the 
long  wavelength  limit,  the  negative  dispersion  from  the  di¬ 
polar  contribution  will  always  dominate  that  from  exchange. 
However,  since  the  prefactor  of  the  logarithm  is  proportional 
to  (kR)2,  for  nanowires  of  small  radii  the  dipolar  term  will 
assert  itself  only  at  very  small  values  of  the  wave  vector.  One 
sees  easily  that  the  dipolar  term  will  dominate  the  long 
wavelength  form  of  the  dispersion  relation  only  when  k 
<(2IR)exp(—DhrMsR2).  Thus,  while  for  m=  1,  we  expect 
the  minimum  frequency  of  the  magnetostatic  mode  to  always 
lie  away  from  zero  wave  vector,  but  in  numerical  calcula¬ 
tions  for  nanowires  of  small  radii,  the  initial  negative  disper¬ 
sion  will  be  a  subtle  feature.  At  larger  values  of  the  radius, 
we  will  see  the  negative  dispersion  clearly. 

For  the  modes  with  m>  2,  from  Eq.  (6b),  we  see  that  the 
dipolar  contributions  to  the  dispersion  provide  negative  dis¬ 
persion  in  the  form  of  a  term  quadratic  in  the  wave  vector, 
with  prefactor  scaling  also  as  R2.  Since  an  exchange  contri¬ 
bution  Dk 2  is  independent  of  radius,  we  expect  to  see  posi¬ 
tive  dispersion  in  these  modes  for  nanowires  of  small  radii, 
and  then  for  nanowires  of  larger  radii,  we  will  realize  nega¬ 
tive  dispersion  initially,  with  the  minimum  frequency  of  a 
given  branch  shifted  away  from  zero  wave  vector. 

One  may  cast  the  discussion  of  the  spin-wave  dispersion 
relation  in  terms  of  two-dimensionless  parameters  h 
=  Hq/4ttM s  and  p  =  l/R,  and,  where  1  =  (D/4ttM  s)112  is  an 
exchange  length.  For  Ni,  the  exchange  length  is  roughly  58 
A,  while  for  Fe  it  is  35  A.  In  what  follows,  we  use  the 
reduced  frequency  to  =  H/4ttMs  ,  and  the  frequencies  will  be 
plotted  against  the  logarithm  of  kR. 

In  Fig.  1(a),  for  m=  1  and  for  h  =  0.19  (this  corresponds 
to  an  applied  dc  field  of  4  kG  for  Fe),  we  show  the  spin  wave 
dispersion  for  the  case  where  p  =  2.  Within  the  range  of  fre¬ 
quencies  displayed,  there  is  one  mode  only,  which  in  the 
limit  of  zero  wave  vector  approaches  the  FMR  frequency, 
which  in  our  dimensionless  units  has  the  frequency  h  + 1/2. 
The  negative  dispersion  at  small  wave  vectors  discussed 
above  does  not  show,  on  the  scale  of  this  plot.  Clearly,  by  the 
time  the  parameter  kR  approaches  unity,  the  excitation  en- 
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(a) 


Log(kR) 

(b) 


Log(kR) 

FIG.  1.  For  the  case  where  m  =  1,  and  for  Ii  =  H0/4ttM s 
=  0.19,  we  show  the  spin- wave  frequencies  plotted  as  a  function  of 
In (kR)  for  (a)  p  =  2  and  (b)  p  =  0.2,  where  p  is  the  ratio  of  the 
exchange  length  defined  in  the  text  and  the  radius  of  the  nanowire. 

ergy  of  the  spin  wave  is  dominated  by  exchange  for  nanow¬ 
ires  of  such  small  radius.  In  Fig.  1(b),  for  the  same  applied 
magnetic  field,  and  again  for  m=  1,  we  show  dispersion 
curves  for  nanowires  of  much  larger  radius,  where  p  =  0.2. 
We  now  see  several  standing  wave  exchange  branches, 
whose  frequency  at  zero  wave  vector  is  above  that  of  the 
FMR  mode.  Also,  at  these  larger  radii,  for  the  mode  that 
approaches  the  FMR  frequency  at  long  wavelengths,  the 
negative  dispersion  at  small  wave  vectors  is  now  evident  in 
the  plot. 

We  find  very  interesting  behavior  at  much  larger  radii,  as 
illustrated  in  Fig.  2.  These  calculations  again  are  for  m  =  1 
and  h  =  0.19,  but  now  we  have  p  =  0.05.  In  the  case  of  Ni, 
this  would  describe  a  nanowire  whose  radius  is  a  bit  under 
1200  A.  We  now  have  a  larger  number  of  more  closely 
spaced  exchange  branches,  as  expected,  but  now  notice  that 
there  are  two  exchange  modes  which  lie  below  the  ferromag¬ 
netic  resonance  frequency,  in  the  limit  of  zero  wave  vector. 
The  magnetostatic  mode  that  approaches  the  FMR  frequency 
as  k  — »  0  shows  large  negative  dispersion  at  small  wave  vec¬ 
tors,  until  exchange  takes  over,  and  in  fact  we  see  that  the 
magnetostatic  mode  hybridizes  with  the  low  frequency  ex¬ 
change  modes,  as  the  wave  vector  increases. 

This  set  of  dispersion  curves,  and  the  comparison  between 
Fig.  2  and  Fig.  1(b)  requires  some  comment.  Let  us  consider 
the  basic  solutions  we  encounter  in  the  limit  the  wave  vector 
k  is  zero.  If  we  think  of  infinitely  extended  ferromagnetic 
matter  for  the  moment,  then  we  shall  have  spin  waves  that 
propagate  within  the  xy  plane,  whose  frequency  is  given  by 
the  well  known  expression  I7(  k)  =  \(Hq  +  D  k2)(Hq 
+  47 tM s+ D k2)]U2,  where  we  denote  them  in  plane  wave 
vector  as  k.  When  such  a  wave  is  confined  to  within  a  cyl¬ 


FIG.  2.  For  m—  1,  fi  =  0.19,  and  p  —  0.05,  we  show  the  spin- 
wave  frequencies  of  the  nanowire  plotted  as  a  function  of  ln(kR). 

inder  of  finite  radius  R,  k  will  obtain  quantized  values  in  the 
range  of  k=  Kn  =  nir/R,  where  the  lowest  standing  mode  has 
n=  1 .  There  is  in  addition  a  uniform  mode  of  the  cylinder, 
wherein  the  magnetization  engages  in  perfectly  circularly  po¬ 
larized  precession  at  the  FMR  frequency  H0+2ttM s  , 
which,  it  will  be  noted,  lies  higher  in  frequency  than  ft(/c 
=  0).  In  the  cylinder,  in  the  presence  of  the  boundaries,  the 
boundary  conditions  mix  the  uniform  FMR  mode  with  the 
standing  wave  exchange  modes,  so  to  speak.  Now  if  the  ra¬ 
dius  of  the  nanowire  is  so  small  that  D, (k  =  tt/R)  lies  above 
the  FMR  frequency,  we  have  a  situation  where  one  of  the 
modes,  which  we  have  referred  to  above  as  the  magnetostatic 
mode,  approaches  the  FMR  frequency  as  k— > 0,  while  all  the 
standing  wave  exchange  modes  lie  above  this  frequency. 
This  is  the  case  for  the  example  given  in  Fig.  1(b).  Now  as 
the  radius  of  the  cylinder  is  increased,  one  or  more  of  the 
modes  of  frequency  IK  k„)  will  drop  below  the  FMR  fre¬ 
quency.  In  the  limit  k— +0,  we  thus  have  standing  wave  ex¬ 
change  modes  which  lie  below  the  FMR  frequency.  This  is 
the  case  for  the  example  in  Fig.  2.  The  magnetostatic  mode, 
which  approaches  the  FMR  frequency  in  this  limit,  exhibits 
negative  dispersion  at  small  wave  vectors  as  expected  from 
Eq.  (6a),  and  the  boundary  conditions  lead  to  the  hybridiza¬ 
tion  with  the  exchange  modes  evident  in  Fig.  2. 

The  situation  just  described  is  very  different  than  realized 
in  the  much  studied  case  of  the  thin  film  magnetized  parallel 
to  its  surfaces.  We  still  have  the  in  plane  standing  wave 
exchange  modes,  with  the  frequency  0,(k)  just  as  above,  but 
now  the  uniform  mode  that  mixes  with  these  is  elliptically 
polarized,  with  frequency  [ //0( //(l  +  4  ttM s)  ] 1  /2=  O ( 0 ) .  The 
standing  wave  spin  waves  always  lie  above  the  ferromag¬ 
netic  resonance  frequency. 

These  comments  have  implications  for  ferromagnetic 
resonance  studies  for  nanowires,  where  the  response  of  the 
sample  is  probed  at  fixed  frequency,  and  the  external  dc 
magnetic  field  is  swept.  In  the  thin  film,  the  standing  wave 
exchange  resonances  are  always  found  at  resonance  fields 
below  the  main  FMR  resonance  field  (provided,  of  course, 
there  are  no  perturbations  at  the  boundary  so  large  as  to 
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FIG.  3.  For  m  =  2,  h  =  0.19  and  (a)  p  =  0.2  and  (b)  p  =  0.05,  we 
show  the  spin  wave  frequencies  plotted  as  a  function  of  In (kR). 

produce  spin  waves  tightly  bound  to  the  surface).  It  will  be 
the  case  as  well  that  in  nanowires  of  small  radius,  the  stand¬ 
ing  wave  resonances  will  lie  lower  in  field  than  the  main 
FMR  line.  However  as  the  radius  of  the  nanowire  is  in¬ 
creased,  one  reaches  a  point  where  one  may  find  standing 
wave  resonances  of  the  structure  at  fields  above  the  FMR 
field;  such  lines  will  be  produced  by  the  exchange  modes 
that  have  dropped  below  the  FMR  frequency. 

In  Fig.  3  and  in  Fig.  4,  we  show  dispersion  curves  for 
modes  with  m  =  2  and  m  =  3,  respectively.  One  sees  that  for 
the  samples  of  small  diameter,  the  dispersion  is  positive  at 
long  wavelengths,  whereas  one  sees  negative  dispersion  at 
larger  radii,  as  expected  from  the  discussion  above. 

B.  The  microwave  response  of  conducting  ferromagnetic 
nanowires 

We  next  turn  to  studies  of  the  microwave  response  of 
ferromagnetic  nanowires,  utilizing  the  formalism  given  in 
Sec.  II B.  All  computations  below  use  parameters  appropriate 
to  Ni  at  room  temperature,  with  the  experiments  of  Ref.  1  in 
mind.  We  take  MS  =  480G,  the  g  factor  to  be  2.15,  the  spin- 
wave  exchange  stiffness  as  D  =  2X  1CT9  G-cm2  and  for  the 
purposes  of  including  the  influence  of  the  skin  depth  on  the 
response,  the  electrical  resistivity  is  7  mi  Ohm  cm.  This  gives 
a  classical  skin  depth  in  the  range  of  a  micron  at  10  GHz, 
which  is  very  much  larger  than  the  radii  of  the  nanowires 
considered  here. 

For  the  purposes  of  calculating  the  FMR  spectra,  we  con¬ 
sider  a  quantity  F  defined  as  follows.  As  we  have  seen  in 
Sec.  II B,  the  precessing  magnetization  of  the  wire  creates  a 
strong  magnetic  field  outside  the  wire,  which  falls  off  in¬ 
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FIG.  4.  Form  =  3,  h  =  0.19,  and  (a)  p  =  0.2  and  (b)  p  =  0.05,  we 
show  the  spin  wave  frequencies  plotted  as  a  function  of  ln(A'R). 

versely  with  the  square  of  the  distance  from  the  axis  of  sym¬ 
metry,  as  we  see  from  Eqs.  (36a)  and  (36b).  The  coupling  of 
the  nanowire  to  a  detector  will  scale  as  the  square  of  the 
strength  of  this  field.  We  define  T  to  be  r  =  ((hp)2 
+  (h<f>)2)P=R+  /h In  where  the  angular  brackets  denote  an  av¬ 
erage  over  the  circumference  of  the  wire. 

Before  we  present  our  results,  we  comment  on  the  role  of 
the  spin  depth  in  the  present  analysis.  Our  point  is  that  it 
differs  qualitatively  for  nanowires  from  the  case  of  the  thin 
film.  First  of  all,  in  a  conducting  magnetic  medium,  the  skin 
depth  is  influenced  not  only  by  the  conductivity  of  the  ma¬ 
terial,  but  by  its  magnetization  as  well.  To  good 
approximation,9  in  place  of  the  classical  skin  depth  <50 
=  cI(  2ttij)(t)  i/2,  penetration  of  microwave  fields  is  con¬ 
trolled  by  SM  =  cl(2TTU)cr/j,v)m,  where  /jlv=(/j.2  — /jl\)I 
is  the  Voigt  susceptibility.  Here  /jl1  =  1  +4irx\ ,  and  fi2 
=  4'7LV'2>  with  X  i  and  Xi  the  dynamic  susceptibilities  defined 
in  Eq.  (2).  For  the  picture  used  here  in  Sec.  II  A,  one  has 
nv=(B20-n2)/(H0B0-n2). 

Now  if  we  wish  to  consider  the  ferromagnetic  resonance 
response  of  thin  films,  note  that  their  resonance  frequency  is 
(HqB o)1/2-  At  precisely  this  frequency,  the  Voigt  susceptibil¬ 
ity  diverges,  if  we  examine  the  simple  expression  above.  A 
consequence  is  the  skin  depth  collapses  to  zero.  Of  course, 
dissipative  effects  not  included  in  the  simple  expression  just 
given  limit  the  divergence,  but  it  is  the  case  that  the  skin 
depth  collapses  to  a  small  value,  as  one  scans  through  reso¬ 
nance.  In  a  high  quality  ferromagnetic  metal  such  as  Fe,  the 
effective  skin  depth  on  resonance  can  be  as  small  as  500  A, 
in  the  10  GHz  frequency  range.  The  very  strong  frequency 
dependence  of  the  skin  depth  as  one  scans  through  resonance 
means  it  is  essential  to  include  its  role,  in  the  discussion  of 
microwave  response  of  thin,  conducting  ferromagnetic  films. 
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In  Ref.  9  and  earlier  references  cited  therein,  one  sees  strong 
influences  of  the  frequency  dependence  of  the  skin  depth, 
including  the  “window”near  B0  where  the  skin  depth  be¬ 
comes  very  large. 

In  the  ferromagnetic  nanowires,  we  are  concerned  with 
much  higher  frequencies,  near  H0  +  2  ttM  s ,  where  there  is 
no  particular  resonance  or  strong  frequency  dependence  in 
the  Voigt  susceptibility.  The  penetration  depth  of  the  fields 
differs  only  nominally  for  that  expected  in  the  absence  of  the 
magnetic  response.  Thus,  for  nano  wires  with  radius  small 
compared  to  the  classical  skin  depth,  their  finite,  metallic 
conductivity  has  at  most  a  modest  influence  on  their  re¬ 
sponse. 

Further  comments  along  these  lines  are  of  interest.  To 
return  to  the  case  of  the  thin  metallic  film,  if  the  film  has 
thickness  small  compared  to  the  renormalized  skin  depth  on 
resonance,  a  situation  encountered  commonly  in  the  study  of 
the  FMR  of  ultrathin  films,  one  might  suppose  that  within  the 
film,  the  microwave  exciting  field  is  uniform  to  excellent 
approximation.  This  is  not  the  case,  if  the  film  is  illuminated 
by  a  microwave  field  incident  from  one  side.  In  such  a  case, 
it  is  an  elementary  matter  to  see  that  within  the  film,  if  its 
thickness  d  is  thin  compared  to  the  skin  depth,  and  also  the 
skin  depth  is  thin  compared  to  a  vacuum  wavelength,  the 
spatial  variation  h{z)  of  the  microwave  magnetic  field  within 
the  films  is  well  approximated  by,  with  k0  =  D,/c, 


h(z)  =  2li 


(1  -z/d)  +  k0S2M/2d 
1  +  k082Mld 


(37) 


Thus,  even  if  the  film  is  very  thin  compared  to  the  skin 
depth,  the  field  drops  linearly  from  a  value  very  close  to  2  h0 
at  the  surface  exposed  the  microwaves,  to  hQkQ82Mld  at  the 
back  surface;  the  field  is  highly  non  uniform  within  the  film, 
until  one  reaches  thicknesses  where  k082M/d  is  of  order  unity. 

The  physical  origin  of  the  highly  nonuniform  magnetic 
field  distribution  within  a  film  illuminated  from  one  side, 
even  though  its  thickness  may  be  small  compared  to  the  skin 
depth,  is  the  poor  impedance  mismatch  between  the  electro¬ 
magnetic  disturbance  within  the  film,  and  the  transmitted 
wave  on  the  output  surface.  The  magnetic  field  thus  drops 
from  a  value  roughly  twice  that  of  the  incident  field  on  the 
illuminated  surface,  to  a  very  small  value  on  the  backside  of 
the  film.  Quite  in  contrast  to  this,  our  conducting  nanowire  is 
embedded  in  a  spatially  uniform  magnetic  field,  so  the  pro¬ 
file  within  the  wire  is  not  influenced  by  such  considerations. 
If  the  wire  radius  is  small  compared  to  the  classical  skin 
depth,  the  exciting  field  the  wire  will  be  rather  uniform  in 
character. 

We  now  present  a  summary  of  our  studies  of  the  micro- 
wave  response.  The  principal  conclusions  can  be  summa¬ 
rized  quite  briefly. 

First  of  all,  we  have  seen  that  the  nanowire  admits  a  rich 
spectrum  of  exchange  dipole  modes  in  the  vicinity  of  k  =  0 
for  the  modes  with  azimuthal  quantum  number  m=  1.  We 
have  modes  both  above  and  below  that  of  the  ferromagnetic 
resonance  frequency,  for  nanowire  diameters  in  the  range  of 
a  few  hundred  Angstroms  or  more.  However,  if  we  assume 
that  the  surface  anisotropy  Ks  vanishes,  then  the  only  mode 


(b) 


FIG.  5.  We  illustrate  the  role  of  surface  anisotropy  in  activating 
the  standing  wave  exchange/dipole  satellite  modes  in  the  ferromag¬ 
netic  resonance  spectrum  of  the  nanowire.  The  calculations  are  for 
a  Ni  nanowire  with  radius  1000  A,  and  focus  on  the  low  field  wing 
of  the  absorption  line.  We  calculate  the  derivative  with  respect  to  dc 
field  of  the  ratio  T  defined  in  the  text  for  (a)  the  case  where  the 
surface  anisotropy  is  zero,  and  (b)  where  Ks=  —0.3  ergs/cm2.  The 
two  insets  show  the  inverse  of  the  determinant  of  the  matrix  that 
must  be  inverted  to  solve  for  the  magnetization  amplitude.  The 
peaks  are  the  positions  of  the  spin  wave  normal  modes.  We  have 
used  the  rather  small  value  G  =  0.5X  108  sec~’  for  the  Gilbert 
damping  function,  to  enhance  the  satellite  features. 


we  see  in  the  calculated  spectmm  is  the  main  ferromagnetic 
resonance  mode  itself.  We  have  searched  even  for  weak 
structures  from  exchange/dipole  modes  removed  from  the 
main  FMR  mode,  to  find  no  evidence  of  absorption  features. 
It  is  very  much  as  if  there  is  a  hidden  theorem  we  have  not 
been  able  to  prove  analytically  which  forbids  all  modes  other 
than  the  FMR  mode  to  be  active  in  microwave  excitation 
when  Ks=  0.  It  is  well  known  that  in  the  elementary  Heisen¬ 
berg  ferromagnet,  the  standing  spin  wave  resonances  are  “si¬ 
lent”  in  the  absence  of  surface  spin  pinning  produced  by 
surface  anisotropy,  in  the  limit  the  microwave  exciting  field 
may  be  viewed  as  uniform  within  the  film.  We  were  sur¬ 
prised  to  see  the  complete  absence  of  all  modes  except  the 
main  FMR  mode  in  our  numerical  studies  of  the  ideal 
nanowire,  with  surface  anisotropy  absent.  In  our  minds,  the 
existence  of  a  selection  rule  is  not  obvious. 

We  illustrate  these  remarks  in  Fig.  5,  which  shows  a  por- 
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tion  of  the  low  field  wing  of  the  main  FMR  line  of  a  Ni 
nano  wire  1000  A  in  diameter.  We  have  assumed  the  fre¬ 
quency  is  34.4  GHz,  as  utilized  in  Ref.  1.  Thus,  the  reso¬ 
nance  field  is  a  bit  above  8  kG,  for  our  parameters.  In  Fig. 
5(a),  we  show  a  region  of  the  low  frequency  wing,  calculated 
for  the  case  where  surface  anisotropy  is  absent.  One  can 
perceive  no  structure  whatsoever  from  exchange/dipole 
modes.  Yet  modes  exist  in  this  field  regime,  and  illustrated  in 
the  inset,  where  we  plot  the  inverse  of  the  determinant  that 
must  be  evaluated  to  find  the  various  field  amplitudes  in  Sec. 
II B.  This  determinant  has  a  peak  at  the  held  where  an  ex¬ 
change  dipole  mode  occurs. 

If  surface  anisotropy  is  introduced,  then  the  exchange  di¬ 
pole  modes  “light  up”  as  illustrated  in  Fig.  5(b).  The  results 
displayed  here  are  for  the  choice  Ks=  —0.3  ergs-cm2,  which 
is  a  substantial  value  of  the  surface  anisotropy.  Spin  pinning 
to  this  degree  may  be  realized  at  the  surface  if,  for  example, 
an  oxide  layer  is  present.  The  principal  satellite  at  the  higher 
held  has  an  intensity  of  about  1%  of  that  of  the  main  FMR 
line. 

In  situations  where  the  nanowire  radius  is  such  that  an 
exchange/dipole  mode  lies  quite  close  in  frequency  to  the 
FMR  mode,  the  mixing  of  the  two  by  the  boundary  condi¬ 
tions  can  endow  the  satellite  with  appreciable  integrated 
strength.  This  is  illustrated  in  spectra  taken  by  Ebels  and 
Wigen.1  These  authors  explore  the  FMR  spectra  of  Ni 
nanowires  with  radius  350,  800,  2700,  and  5000  A.  With  the 
exception  of  the  800  A  sample,  a  single  absorption  line  is 
observed.  However,  for  the  800  A  sample,  there  is  a  clear 
doublet.  In  Fig.  6,  we  show  calculations  of  the  spectrum  for 
three  of  the  sample  radii,  and  indeed  we  find  a  doublet  very 
similar  in  character  to  that  observed.  The  relative  oscillator 
strength  of  the  two  modes  is  given  nicely,  when  theory  is 
compared  to  experiment.  Not  displayed  is  the  calculation  for 
the  5000  A  case,  which  also  shows  a  single  feature  as  in  the 
data.  The  calculations  in  Fig.  6  use  Ks=  —0.85  ergs/cm2, 
corresponding  to  a  hard  axis  normal  to  the  nanowire  surface. 
This  is  again  a  large  value  of  the  surface  anisotropy,  appro¬ 
priate  to  a  case  where  oxide  is  present  on  the  surface. 


IV.  CONCLUDING  REMARKS 

We  have  developed  the  theory  of  the  exchange/dipole 
spin  wave  modes  of  ferromagnetic  nanowires  of  circular 
cross  section,  along  with  that  of  their  response  to  spatially 
uniform  microwave  exciting  fields.  We  also  have  presented 
studies  of  both  the  excitations  in  and  the  microwave  response 
of  samples  quite  similar  to  those  explored  in  Ref.  1,  and 
compared  the  differences  between  the  wire  and  the  well 
know  case  of  the  thin  film.  We  conclude  with  some  addi¬ 
tional  remarks. 

The  origin  of  the  linewidth  in  such  samples  is  of  interest. 
Of  course,  there  is  the  intrinsic  linewidth,  described  in  our 
phenomenology  by  the  damping  term  in  the  Landau  Lifshitz 
equation  of  motion.  In  addition,  there  will  be  extrinsic 
mechanisms  operative  as  well.  For  the  case  of  ultrathin  fer¬ 
romagnetic  films,  in  recent  work  we  have  shown6  that  the 
two  magnon  mechanism  considered  many  years  ago  as  the 
source  of  extrinsic  linewidth  in  the  garnets9  can  operate  in 
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FIG.  6.  For  three  values  of  the  radius  of  Ni  nanowires,  (a)  350  A 
(b)  800  A,  and  (c)  2700  A,  we  show  calculated  ferromagnetic  reso¬ 
nance  spectra.  The  calculations  employ  G  =  2.5X  108  sec-1  which  is 
a  realistic  value  for  Ni,  and  we  have  taken  Ks=  —0.85  ergs/cm2. 


the  ultrathin  films  as  well,  by  virtue  of  the  fact  that  in  two 
dimensions,  the  long  ranged  dipolar  interaction  in  combina¬ 
tion  with  exchange  produces  an  off  center  minimum  in  the 
spin  wave  dispersion  relation,  for  a  range  of  propagation 
directions.  This  leads  to  spin  wave  modes  of  finite,  nonzero 
wave  vector  degenerate  with  the  zero  wave  vector  FMR 
mode.  Surface  defects  (or  defects  of  any  other  sort)  may  then 
scatter  energy  from  the  FMR  mode  to  the  finite  wave  vector 
modes,  producing  a  dephasing  contribution  to  the  linewidth. 
Our  predictions  seem  confirmed  by  recent  experiments.7  One 
can  inquire  if  the  two  magnon  mechanism  is  operative  in 
nanowires  as  well. 

In  principle,  this  mechanism  should  be  operative.  First 
suppose  we  have  a  nanowire  with  undulating  surfaces  of 
such  character  that  the  wire  is  still  form  invariant  under  ro¬ 
tations  about  the  z  axis.  Then  from  symmetry  considerations, 
there  will  be  a  nonzero  matrix  element  for  scattering  energy 
from  the  zero  wave  vector  FMR  mode  of  m  =  1  character,  to 
m  =  1  modes  of  finite  wave  vector  on  the  FMR  branch  which 
are  degenerate  with  the  FMR  mode.  We  see  when  Eq.  (6a)  is 
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supplemented  by  an  exchange  term  that  scales  as  Dk1 2 3 4 5,  there 
is  always  a  minimum  in  the  dispersion  relation  of  m  =  1  spin 
waves  away  from  zero  wave  vector,  and  thus  there  always 
will  be  finite  wave  vector  modes  degenerate  with  the  FMR 
mode.  However,  when  R<RC=(D/ ttM s)l/2,  from  the  struc¬ 
ture  of  the  dispersion  curve  the  minimum  will  lie  quite  near 
zero  wave  vector,  and  the  number  of  degenerate  modes  will 
be  rather  small.  When  R>RC,  the  minimum  will  be  in  the 
vicinity  of  k=  HR,  and  we  expect  the  mechanism  to  be  more 
efficient.  For  Ni,  Rc  is  roughly  100  A,  and  for  Fe  perhaps  a 
factor  of  2  shorter. 

If  the  surface  has,  say,  roughness  of  random  character,  so 
that  azimuthal  symmetry  is  no  longer  present,  then  the  ma¬ 
trix  element  for  scattering  energy  from  the  m  =  1  FMR  mode 
to  modes  with  in >  1  is  nonzero.  The  criterion  R>RC  re¬ 
mains  relevant,  in  that  it  is  for  these  larger  radii  that  the 
dispersion  relation  for  the  m  >  1  magnetostatic  modes  (those 
which  approach  the  FMR  frequency  as  the  wave  vector  van¬ 
ishes)  has  an  off  center  minimum.  Also,  this  is  a  crude  cri¬ 
terion  for  the  appearance  of  zero  wave  vector  modes  below 
the  FMR  frequency  of  the  cylinder,  which  is  another  source 
of  finite  wave  vector  modes  degenerate  with  the  FMR  mode. 

The  above  discussion  suggests  that  the  two  magnon 
mechanism  should  be  relatively  inefficient  in  very  small  ra¬ 
dius  nano  wires  with  R<RC,  but  for  larger  nanowires  we 
suggest  it  should  be  quite  efficient.  A  detailed  theory,  beyond 
the  scope  of  the  present  paper,  will  be  required  to  explore 
this  question  further. 

Another  aspect  of  the  geometry  considered  here  is  the 
large  microwave  magnetic  field  created  outside  the  wire,  by 
the  precession  of  the  magnetization.  We  see  from  Eq.  (18) 
that  for  a  spin  wave  with  finite  wave  k,  we  have  a  field  that 
falls  off  as  exp(  —  Kp)/pV2  far  from  the  wire.  For  long  wave¬ 
length  modes,  this  field  has  a  very  long  range.  In  the  limit  of 
zero  wave  vector,  we  see  from  Eqs.  (36)  that  the  field  created 
by  the  precession  of  the  magnetization  falls  off  inversely 
with  the  square  of  the  distance  from  the  center  of  the  wire. 
The  existence  of  this  large  field  outside  the  nanowire  is  a 
substantial  difference  that  the  situation  with  thin  films.  If,  for 


the  thin  film,  we  consider  the  uniform  precession  mode  ex¬ 
cited  in  an  FMR  experiment,  the  macroscopic  magnetic  field 
created  by  precession  of  the  magnetization  is  completely 
confined  to  within  the  film.  There  are  “magnetic  poles”  on 
the  film  surfaces  associated  with  the  precessing  magnetiza¬ 
tion,  but  the  field  generated  by  these  is  confined  to  within  the 
film  itself,  as  in  a  condensor  plate  problem.  If  we  consider  a 
spin  wave  whose  wave  vector  parallel  to  the  film  surfaces  is 
k,  there  is  a  macroscopic  field  outside  with  the  spatial  varia¬ 
tion  exp(— kz)  if  the  film  surfaces  are  parallel  to  the  xy  plane, 
but  the  prefactor  which  controls  the  strength  of  this  field 
scales  as  AttM  s{kd),  where  d  is  the  film  thickness.  This 
field  is  thus  very  weak  in  the  thin  film  limit,  or  whenever  the 
wavelength  of  the  spin  wave  is  large  compared  to  the  film 
thickness. 

The  existence  of  these  large  fields  outside  the  nanowire 
have  interesting  implication  for  samples  such  as  those  ex¬ 
plored  by  the  authors  of  Ref.  1.  We  have  long  ranged  inter¬ 
actions  between  the  nanowires,  so  one  is  led  to  explore  the 
collective  excitations  of  the  nanowire  array.  In  the  particular 
case  of  the  samples  used  by  Ebels  and  Wigen,  the  nanowires 
are  accurately  parallel  to  each  other,  but  randomly  arrayed 
over  a  plane.  We  thus  have  a  random  two-dimensional  (2D) 
system  with  long  ranged  interactions.  Indeed,  if  one  were  to 
excite  spin  waves  of  variable  wavelength  k,  as  in  a  light 
scattering  experiment,  then  the  ratio  of  the  average  interwire 
separation  to  the  range  of  interaction  can  be  varied.  There  is 
an  analogy  between  this  system,  and  the  vortex  glass  state  of 
2D  superconductors,  where  here  the  average  vortex  separa¬ 
tion  can  be  varied,  but  the  range  of  the  interaction  if  fixed  as 
the  London  penetration  depth.  In  the  case  of  nanowire  arrays, 
the  nature  of  their  collective  excitations  will  be  a  most  inter¬ 
esting  topic,  for  ordered  and  disordered  arrays. 
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